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Using a novel approach, we investigate the shape of the av- 
erage spectrum and the spectral fluctuations of the fc-body 
embedded unitary ensemble in the limit of large matrix di- 
mension. We identify the transition point between semicircle 
and Gaussian shape. The transition also affects the spectral 
fluctuations which deviate from Wigner-Dyson form and be- 
come Poissonian in the limit k -C m -C I. Here m is the 
number of Fermions and I the number of degenerate single- 
particle states. 

PACS numbers: 05.30.-d, 21.60.Cs, 31.10.+Z, 73.20.Dx 



Introduction. The stochastic behavior displayed by 
spectra and wave functions of quantum many-body sys- 
tems (atoms, molecules, atomic nuclei, quantum dots) is 
usually and successfully modelled in terms of canonical 
random-matrix theory (RMT) However, this type 

of modelling is not completely realistic: All the above- 
mentioned many-body systems are effectively governed 
by one- and two-body forces, while canonical RMT is 
tantamount to assuming many-body forces between the 
constituents. Thus, a stochastic modelling of the one- 
and two-body interaction would yield a much smaller 
number of independent random variables than used in 
RMT. For instance, the number of independent two- 
body matrix elements in a shell-model calculation in 
atoms or nuclei is typically much smaller than the di- 
mension N of the matrices involved, while the number 
of independent random variables in RMT is of order N 2 . 
This difference poses the question whether a more real- 
istic stochastic modelling of many-body systems might 
yield results which differ from RMT predictions. The 
question was addressed in the 1970's with the help of nu- 
merical simulations using matrices of fairly small dimen- 
sions. The main results were: In a certain limit, the aver- 
age level density does not have the shape of a semicircle 
but is Gaussian; the ensembles are neither stationary nor 
ergodic; unfolding of the spectra yields Wigner-Dyson 
spectral fluctuation properties, see Refs. @. Interest in 
model Hamiltonians with random two-body interactions 
has resurged in recent years in several areas of many- 
body physics (see Ref. [Q and further references therein), 
and the question of possible further differences between 
such models and RMT has resurfaced. It is the purpose 
of this letter to shed new light on this question. In par- 
ticular, we will show that for certain parameter values, 



the spectral fluctuation properties of realistic models dif- 
fer from those of RMT. We discuss the implications of 
this novel and surprising result. 

Focussing attention on the case of unitary symme- 
try, we use a paradigmatic model, the fc-body embed- 
ded Gaussian random ensemble EGUE(fc) introduced by 
Mon and French ||. For this ensemble, a number of an- 
alytical results has been obtained [||,fi]||. However, as 
emphasized in Ref. Q, page 413, a comprehensive an- 
alytical approach to the spectral fluctuation properties 
of the ensemble is still lacking. This is due to the fact 
that, in contrast to RMT, embedded ensembles do not 
possess the (orthogonal or unitary) invariance in Hilbcrt 
space which is so essential for the successful analytical 
treatment of RMT. 

We use three different methods: The supersymmetry 
approach, the "binary correlation approximation" of Mon 
and French || , and the construction of two "limiting en- 
sembles". We also make use of two fundamental and 
novel results on the second moment: The eigenvector ex- 
pansion and duality. We compare our results with those 
of the Gaussian unitary ensemble of random matrices 
(GUE). 

Definitions. Our Hilbert space is spanned by N = ( ) 
Slater determinants // = l,...,N, obtained by dis- 
tributing m spinless Fermions over I degenerate single- 
particle states. The ratio / = m/l is the filling factor. 
Using standard creation and annihilation operators Oj 

and a,j, j = 1, . . . , I, we have = Y\™ =1 a] s |0). Here |0) 
is the vacuum state. The creation operators are written 
in ascending order of the indices j s . The fc-body inter- 
action Vk has the form {k = 1, . . . , m) 
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Here and in Eqs. (Q), (|^), the prime on the summation 
sign indicates that the sums over all indices run from 1 to 
I with the constraint that j\ < . . . < jk and i\ < . . . < ik- 
The fc-body interaction matrix is complex 

Hermitean. The independent matrix elements are uncor- 
related Gaussian random variables with zero mean and 
a common second moment v 2 . The value of Vv 2 deter- 
mines the overall energy scale and is set equal to unity 
without loss of generality. By taking matrix elements 
(^|Vfc|/i), we "embed" the random fc-body interaction 
into an m-Fermion system. This defines EGUE(fc). The 
ensemble is invariant under unitary transformations of 
the single-particle states and, for fc = m, reduces 
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to the standard GUE. As always in RMT, we are inter- 
ested in the limit N — > oo (I — > oo). 

Second moment. The Gaussian distribution of the ma- 
trix elements in Eq. (|l|) implies that the (z/|Vfe|/i)'s are 
also Gaussian distributed random variables with zero 
mean value. All spectral properties of EGUE(fc) are, 
therefore, determined by the second moment 
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The bar denotes the ensemble average. We have A$, t - 
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^pa.pv dllU \l\pv.pa) — ^ap,vp — ^pv.pa- -L lie ma- 
trix is Hermitean in the pairs of indices (/i, v) and 
(p, a). Moreover, it is easy to prove the "duality" relation 

A^pa- — -^pa.Jv which connects the second moments of 
the fc-body and the (m — fc)-body interaction. 

Eigenvector expansion. We construct the eigenvectors 
(j( sa ) anc i eigenvalues A^(k) of the Hermitean matrix 
A^ satisfying £ pCT A^df = A«(*)C& o) . Here 
s = 0, . . . , m and a labels the degenerate eigenvectors. 
We are guided by the example of the GUE where the 
second moment of the Hamiltonian H reads Hp a H p[J = 
(A 2 /N)8p U 8 pa . The two Kronecker symbols display the 
unitary invariance of the GUE. The matrix b ap is eigen- 
function of the second moment with eigenvalue A 2 . All 
traceless unitary matrices are likewise eigenfunctions but 
belong to eigenvalue zero. In the present case, we use 

the ansatz cjfu — (t l \ a j 1 ■ • • a ] 3 a i s ■ ■ • a h W)- The label 
a enumerates all possible distinct choices of the indices 
ji < . . . < j s ; i\ < . . . < i s . It is easy to check that 
is eigenvector of the matrix A^ if no two indices (j r , i r t) 
are equal. The corresponding eigenvalue is 
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When at least two indices (j r ,i r >) are equal, cjfu^ is 
not an eigenfunction of A^ but a linear combination 
of eigenfunctions with labels s' < s. This is because 
an eigenfunction Cjf v a ^ with eigenvalue A^ s ^(k), s' < s, 
looks like a member of class s when the operator defin- 
ing C^ s ' a) is multiplied by (£V a p a p ) s ~ s ' = m s ~ s ' . To 
remove the components belonging to lower s'-values, we 
orthogonalize (in the sense of the trace) all matrices C^t 
in which at least two indices (j r , i r >) are equal, to all ma- 
trices generated from classes s' < s in the manner just 
described. The resulting matrices are eigenvectors with 
eigenvalue A^(k). We choose Hermitean linear combi- 
nations of the degenerate eigenvectors which obey the or- 
thonormality condition ^„ v CpfiCvif} — NS st S a b. The 
number of linear independent eigenvectors in class 
s is given by = 1 and 
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We have YlT=o ^ = N 2 showing that the eigenvectors 
form a complete orthonormal set. Hence, the matrix A^> 
possesses an eigenvalue decomposition of the form 



A W = 1 V A< s 

s=0 
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We note that the eigenvectors do not depend on the rank 
k of the interaction, only the eigenvalues do. Eq. (|^) 
shows that the sum over s actually terminates at s = 
m — k. For k = m only s — contributes, and the 
result reduces to the GUE expression with A 2 = A^°\m). 
Conversely, Eq. (§) extends the GUE result to EGUE(/c) 
and constitutes a central result of this paper. 

Moments of Vk ■ Using the eigenvector decomposition 
of A^ k \ duality, and the orthonormality of the C^ a ^'s, 
we calculate the low moments of Vk and the kurtosis k 
for which we write n = 2 + Q(k,m,l). We recall that 
k= 2(3) for the semicircle (Gaussian, respectively). We 

find Q(k,m,l) = (l/JV)Erio (TO_fc ' fc) [( A(s) WA (s) (m - 
k))/(A^{k)) 2 )}D (s K For I -► oo and keeping both k 
and m fixed, we have Q(k,m,l) — ► if 2k > m while 

Q(k,m,l) -> ( m k k )/Ck) for 2fc ^ m - This sll0ws that 
the transition of the average spectrum from semicircle to 
Gaussian shape begins (with decreasing k) at 2k = m. 
We ascribe the special role of the point 2k = m to duality. 
Likewise we can show that the relative fluctuations of the 
first and second moments of Vk and, thus, non-ergodic 
features vanish for I — > oo. 

Super symmetry. We calculate spectrum shape and 
spectral fluctuation properties of EGUE(/c) using the su- 
persymmetry approach After averaging over the 

ensemble, the integrand of the generating functional con- 
tains an exponential whose argument depends linearly 
on the matrix A( k \ With the eigenvalue decomposi- 
tion (Eh, the argument of the exponential becomes a sum 
of squares of bilinear forms in the integration variables. 
This allows us to perform the Hubbard-Stratonovich 
transformation. For each pair (s, a) we introduce a su- 
permatrix er( sa ) of composite variables. The resulting 
integral over the composite variables contains the factor 
exp(— £ c fi). The effective Lagrangean £ e s is given by 
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Here, A (s) (fc) = +[A (s) (fc)] 1/2 . The energy arguments Ei 
and E2 of the advanced and the retarded Green function 
define E = (1/2)(E 1 + E 2 ) and e = E 2 - E u while J 
stands for the source terms. The diagonal supermatrix 
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L distinguishes the retarded and advanced cases and is 
defined in Ref. ||. The saddle-point equation has the 
solution CTg S p.^ = S s qt^ where t(°) is the standard GUE 
saddle-point solution. Thus, the saddle-point condition 
yields a semicircular spectrum and universal GUE spec- 
tral fluctuations. To determine the range of validity of 
this solution, we have calculated the first non-vanishing 
term in the loop expansion. The expansion is obtained 
by writing <j( s °) = Us.p: + Sa^ sa ^ and expanding in pow- 
ers of Sa^ ^ . We recall that for the GUE, each term of 
the loop expansion vanishes in the limit N — > oo with 
an inverse power of N. For the one-point function, we 
find that the loop correction is proportional to Q(k, m, I). 
This is consistent with the result of the previous para- 
graph and reaffirms our conclusion that the transition 
from semicircular to Gaussian shape sets in at 2k = m. 
It would be desirable to show that all higher terms of the 
loop expansion vanish likewise asymptotically for 2k > m 
but this proof is beyond our means. For the two-point 
function, the loop correction yields non-universal spec- 
tral fluctuations of the type first considered by Kravtsov 
and Mirlin The amplitude of this correction (given 
in units of the inverse mean level spacing) vanishes, how- 
ever, for I — > co. It does so like N~~ 2 for k = m but only 
like (h\N)~ 2k for both k and / fixed. We see that for 
I —t oo, the supersymmetry approach does not yield a 
limit on the range of validity of the universal GUE spec- 
tral fluctuations of EGUE(fc). 

The case k <C m -C I. For this case where the su- 
persymmetry method does not yield relevant informa- 
tion on spectral fluctuations, we use a modification of 
the "binary correlation approximation" In the av- 

erage two-point function g(zi)g(z2), we expand both 
traced Green functions g in powers of V& . We collect the 
terms containing equal powers of Vfc. The ensemble av- 
erage is taken by Wick-contracting pairs of Vfc's in all 
possible ways. We evaluate all pairs located on the same 
Green function as in Ref. ||. For the rest, we use simple 
counting arguments. We show that for k -C rn <C /, terms 
where s pairs of Vk 's are not on the same Green function, 
are smaller by at least a factor l~ sk than the terms with 
s = 0. Thus, the connected part of g(zi)g(z2) vanishes 
asymptotically compared to the disconnected part, and 
the two-point correlation function 
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approaches zero in the limit k <C m <C I. For similar rea- 
sons, all higher correlation functions also vanish in the 
same limit, and the spectral fluctuations become Poisso- 
nian. 

The case k = 1 with k <C to <C Z illustrates this result. 
Taking first m = k = 1, we construct the eigenvalues e a 
and eigenfunctions \I/ Q of a given realization of EGUE(l) 
by diagonalization. The e Q 's obey GUE spectral statis- 



tics, the average spectrum has semicircle shape and the 
coefficients U a j in the expansion = 53j=i ^aj a \\^) 
in terms of the single-particle basis are Gaussian dis- 
tributed random variables. For to > 1 and every re- 
alization of the ensemble, the eigenfunctions are Slater 
determinants x °f the ^ a 's, the eigenvalues are sums of 
the e a 's, and the average spectrum is an m-fold convo- 
lution of the semicircle. This shows immediately that for 
m>l, the spectrum has Gaussian shape, and that the 
spectral fluctuations are Poissonian. We have also shown 
that the states x are localized. 

Regular graphs. Further insight into the spectral prop- 
erties of EGUE(fc) is gained by using yet another ap- 
proach involving regular graphs and limiting ensembles. 
A graphical representation of EGUE(fc) is obtained by 
assigning to each Hilbert-space vector a vertex fi, 
and to each non-diagonal element (z/|T4|/i) which does 
not vanish identically, a link connecting the vertices v 
and fi. The diagonal matrix elements (/i|T4|/j) are repre- 
sented by loops attached to the vertices fi. The number 
of vertices is N. The number M of links emanating from 
a given vertex is the same for all vertices and given by 

m =£(>;)('-;»). 

For k < m, we have M < N- 1 while M = N - 1 for k = 
to. The resulting graphical structure is called a "regular 
graph" in the mathematical literature. The total number 
P of links is given by P = (1/2)MN. The number K of 
uncorrelated matrix elements of EGUE(fc) is given by 
K = (f) . It is interesting to study the ratio K/P. For 
fixed rn and / < 1/2, K/P grows monotonically with k, 
starting out with very small values and reaching the limit 
2N/(N — 1) for k = m. This suggests that deviations of 
EGUE(fc) from universal GUE behavior are caused by the 
fact that the number of independent random variables is 
too small in comparison with the number of links, making 
it impossible for the system to become thoroughly mixed. 

Limiting ensembles. To test this hypothesis, we 
have constructed and analyzed two limiting ensem- 
bles. The first, EGUE m ; n (fc), is given by the ma- 
trix elements (H^fc"'"^) °f the interaction V k min = 
vY2' a j 1 ■ ■ ■ a j k a ik ■ ■ ■ a ii- The factor v is a Gaussian com- 
plex random variable. The ensemble EGUE m i„(/c) has 
the same graphical representation as EGUE(/c) but pos- 
sesses only a single random variable. We succeeded to 
prove that for k = 1 and k — to, EGUE m i n (fc) is fully in- 
tegrable, has a Gaussian average spectrum and spectral 
fluctuations which are not of GUE type. We are con- 
vinced that these properties hold for all k. The second 
limiting ensemble, EGUE max (fc), carries the maximum 
number of uncorrelated random variables consistent with 
the graph structure of EGUE(/c) and has Hilbert-space 
matrix elements given by 
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The matrix v v ^ is complex Hermitean. Elements not 
connected by symmetry are uncorrelated complex Gaus- 
sian random variables with mean value zero and vari- 
ance v vtl y v i ^ = 8^1 5 vll i . Using supersymmetry, we have 
shown that EGUE max (fe) has an average spectrum of 
semicircle shape, that the spectral fluctuations are of uni- 
versal GUE type, and that for I — > oo the leading term 
of the loop expansion vanishes both for the one-point 
and for the two-point functions. We conclude that the 
spectral properties of EGUE max (fc) coincide with those 
of GUE. 

Conclusions. We have studied the shape of the av- 
erage spectrum and the eigenvalue fluctuations of the 
embedded ensemble EGUE(fc) in the limit of infinite 
matrix dimension, attained by letting the number I 
of degenerate single-particle states go to infinity. We 
have shown that for sufficiently high rank k of the ran- 
dom interaction (2k > m where m is the number of 
Fermions), EGUE(fc) behaves generically: The spectrum 
has semicircle shape, and the eigenvalue fluctuations 
obey Wigner-Dyson statistics. This does not come as 
a surprise. A transition to a different regime takes place 
at or near 2k = m. It has long been known that the av- 
erage spectrum changes into Gaussian shape, although 
the point of departure from the semicircle shape was 
not known previously. We have presented conclusive ev- 
idence that in addition — and contrary to general ex- 
pectations — the level fluctuations also change and are 
not of Wigner-Dyson type for 2k < m. In the extreme 
case k <C m <C I, the spectral fluctuations are Poisso- 
nian, and the eigenfunctions are likely to display local- 
ization in Fock space. We cannot pin down precisely the 
fc-value where such change occurs nor can we penetrate 
deeply into the intermediate regime. This is not sur- 
prising as we do not know of any other case where such 
an aim would have been achieved. We cannot even say 
definitively whether the transition from Wigner-Dyson 
to Poissonian statistics is smooth or sudden. But we 
have circumstantial evidence for a smooth transition: (i) 
The non-universal fluctuations calculated from the loop 
correction set in smoothly, (ii) In the case k = 1, the 
transition from GUE (for m = 1) to Poisson behavior 
(for k <C m) is smooth, (hi) The ratio K/P of the num- 
ber of uncorrelated random variables over the number of 
links changes smoothly with k for fixed m. 

Formally and using diagrammatic language, we as- 
cribe the deviations from universal GUE behavior to the 
fact that with decreasing k, intersecting Wick contrac- 
tion lines gain increased weight. Universal GUE results 
are obtained whenever such contributions are negligible. 
This apparently is the case for 2k > m. We ascribe the 
special role of the transition point 2k = m to duality. 

Physically, our results can be understood in terms of 
the ratio K/P of the number of uncorrelated random 



variables over the number of links. We have shown that if 
all links were to carry uncorrelated random variables, the 
ensemble would have GUE spectral fluctuations. Con- 
versely, if all links were to carry the same random vari- 
able, the ensemble would be completely integrable and 
display Poissonian statistics. These statements hold for 
all values of k. The actual situation is located between 
these two limits. EGUE(l) is closest to the integrable 
case, and EGUE(m) corresponds to the GUE. This shows 
that deviations from GUE behavior are not caused by the 
number of zeros in the matrix representation of the in- 
teraction but are strictly due to correlations between the 
matrix elements carried by the links. 

As mentioned in the introduction, numerical simula- 
tions have shown good agreement between the spectral 
fluctuations of EGUE(fc) and those of GUE. We ascribe 
this result to the fact that the dimensions of the matri- 
ces used were quite small. It is easy to see that the ratio 
K/P is relatively close to unity for typical values like 
k = 2, rn = 8, and I = 20. It is only in the limit I ^> 1 
that K/P becomes very small, resulting in significant de- 
viations of EGUE(fc) spectral fluctuation behavior from 
Wigner-Dyson statistics. 

Details, further results and the extension to EGOE(fc), 
the orthogonal case, are given in Ref. [{Tcfl , 
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